s-s*-d-wave superconductor on a square lattice and its BCS phase diagram 
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We study an extended Hubbard model with on-site repulsion and nearest neighbors attraction which tries 
to mimic some of the experimental features of doped cuprates in the superconducting state. We draw and 
discuss the phase diagram as a function of the effective interactions among electrons for a wide range of 
doping concentrations. We locate the region which is relevant for the cuprates setting some constraints on 
the parameters which may be used in this kind of effective models. We also study the effects of temperature 
and orthorrombicity on the symmetry and magnitude of the gap function, and map the model onto a simpler 
linearized Hamiltonian, which produces similar phase diagrams. 
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Accurate results using a variety of_gxperimental 
techniques, like, for instance- ARPES,Bd Josephson 
tunneling,cl penetration depthja or thermal conductivity 
measurements^ provide evidence on the d-wave symme- 
try of both the gap and order parameter functions of 
optimally doped and underdoped cuprates. Penetration 
depth measurements in slightly underdoped samplescl de- 
termined that their critical behavior fall on the 3D XY 
universality class (see also Refs. @||). This fact, which 
is not consistent with BCS weak coupling theory, indi- 
cates that the phase transition corresponds to the Bose- 
Einstein condensation of a single, complex order parame- 
ter field. The existence of a pseudogap with d-wave sym- 
metry ahpi£c the superconducting state for underdoped 
cupratesu E3 provides further confirmation on the non- 
BCS nature of the superconducting state. 

There are several theoretical schemes which seem to fit 
into this experimental staterof affairs: (a) In the magnetic 
scenario for the cuprates,t3 ! til the strong on-site repulsion 
among electrons gives rise to antiferromagnetic collec- 
tive excitations. These degrees of freedom dress the bare 
interaction among the residual electrons, providing the 
pairing mechanism for the d-wave superconducting state 
and giving rise to the pseudogap in the normal state. 

(b) In some spin-charge separation theoriesjij the pseu- 
dogap is related to the pairing of spinons and the su- 
perconducting state comes about when holons condense. 

(c) A further line of thought supposes that an undeter- 
mined high-energy pairing mechanism (which might be 
dressed vertex of case (a)) gives rise to an effective low 
energy pairing Hamiltonian where the coupling constants 
set the d-wave superconducting state in an intermediate 
coupling region]!!] 

A thorough study of the phase diagram and qualita- 
tive features of a simple model which retain the electronic 
structure found in ARPES experiments is therefore an 
urgent task, if superconductivity in the cuprates has any- 
thing to do with scenarios (a) or (c). A BCS treatment of 
such a model should serve as a starting point for more ac- 
curate solutions which deal properly with the strong cor- 



relations of the intermediate coupling regime. It should 
also be a useful tool for studies of transport, magnetic or 
optic properties of the cuprates, where bulk or surface 
impurities, and other inhomogeneities must be taken into 
account Ji3ii3 

We present in this article a detailed study of an effec- 
tive pairing model, where the kinetic energy term comes 
from a tight-binding fit to the_-ARPES band structure 
of Bi 2 Sr2CaCu208+,5 (Bi22I2)j23 and we use either an 
attractive or repulsive on-site interaction plus a nearest- 
neighbors attraction term. In other words: we view the 
ARPES band structure of Ref. [2(] as the low-energy pair- 
ing action of the actual Hamiltonian where the strong in- 
teractions among electrons, charge fluctuations, etc, have 
renormalized the kinetic energy and interaction terms 
giving rise to a set of hopping integrals, an effective on- 
site repulsion and a nearest-neighbors attraction, while 
other possible terms are irrelevant. The only channel 
which has not been dealt with in this kind of experi- 
mental implementation of the Renormalization Group is 
the superconducting one, so we need to perform only a 
conventional BCS decoupling scheme (the Hartree term 
would double-count interaction effects already taken into 
account). This model generates pure s and d, as well 
as mixed s+s*, s+s*+d and s+s*+id superconducting 
states in different regions of the phase diagram, where s* 
denotes an extended s state. We have performed most of 
our calculations at zero temperature where it is believed 
that the BCS approximation provides the gross features 
of the true ground sta.to.of the system, at least for small 
enough interactions.cjO 

We assume a tight-binding model on a squared lattice 
with a dispersion relation which keeps hopping terms up 
to five nearest neighbors 

5 

ak) = ^2t im (k)-n, (i) 

i=l 

where the hopping integrals U and functions rji(k) are 
listed in table I of Ref. E0[ The resulting electronic struc- 
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ture does not possess electron-hole symmetry and has a 
band width of about 1.1 eV. We vary /j, to change the 
occupation number n = 1 — x where x is the hole doping 
in the sample. The experimental evolution of the Fermi 
Surface (FS) with doping is not inconsistent with such a 
rigid shift of the band structure even though the ascrip- 
tion of a FSj-tp a system with a pseudogap is a matter of 
controversy.EiJ We have also performed the calculations 
keeping only t\ and t 2 , and obtained results quantita- 
tively similar to those explained in this article. 

The interaction term has an on-site potential, which 
can be either positive or negative plus a nearest neighbors 
attraction 



V mt {kS,) = -V - V 2 (cos(fc x £) + cosfoO). 



(2) 



where £ is some typical length scale. This potential 
resembles the effective pairing interaction found in the 
magnetic scenario of the cuprates, if we set £ equal to 
the magnetic coherence length (and slightly larger than 
the lattice spacing, see Figure 11 of Ref. |l3|) . 

The BCS expressions for the thermodynamic potential 
and gap functions are 

0(T, M )=MlM + M i^l 2 + |A d | 2 , 
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A(k)= A s + A a *F s *(k)+A d F d (k), 



(3) 



where M is the number of sites, Ek — vffe + l^fcl 2 * s 
the dispersion relation for the quasiparticles, and F s * td — 
cos(k x £) ± cos(fc y ^). 

The gap is a function of the dimensionless variable fc£ 
for this particular choice of the interaction term, but we 
believe that such a radial dependence also holds for any 
potential of the form V(fe£), even for moderate values of 
the coupling constants. In the magnetic scenario of the 
cupratesJi-3 £ is the correlation length and, hence, the fur- 
nished gap function is not commensurate with the Bril- 
louin zone and varies along with x. Plots of A(k) in the 
TM and MY directions for sMiie underdoped and opti- 
mally doped Bi2212 samples Jffi if possible, might demon- 
strate if this conjecture is correct. We set nevertheless £ 
equal to the lattice spacing a throughout this paper. 

The saddle point equations can now be easily written 
because the pairing interactions Vo t2 (k— k') are separable: 



A, 



M 



V^tanh(^) 



V 2 y, A I, F. 



2Ei 



^tanh(^) 



(4) 



We compare the free energy per site of the different 
superconducting states f(T, n) = fi(T, fj,(n))/M +fJ.(n) n, 
to ascertain their relative stability. 
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FIG. 1. Phase diagram for (a) x = 0.05 (b) x = 0.2 and 
(c) x — 0.5. Labels I, II and V denote the regions of metallic, 
d-wave and s+s* states, while in regions III and IV both d 
and s+s* solutions are found. The boundaries of the s+s*+id 
state are drawn only in figure (b), with dashed lines. 



We examine first the phase diagram [Vo,Vy for x in 
the range [0.0-0.5] in order to locate the region of pa- 
rameter space which can be used to best represent the 
superconducting state of the cuprates. Figure 1 shows 
the phase diagrams for x=0.05, 0.2 and 0.5. The coupled 
BCS equations possess only the trivial normal state so- 
lution in region I; only a pure d state solution in region 
II; only a mixed s+s* solution in region V, and only a 
pure s solution on the positive Vq axis. There exist both 
pure d and mixed s+s* solutions in regions III and IV; 
the d state being more stable in region III and the mixed 
s+s*, in region IV. We also obtain a mixed s+s*+d solu- 
tion for rather large values of V 2 , and a mixed s+s*+id 
state. These two states are always the most stable as 
long as they exist. The frontiers of stability of the dif- 
ferent solutions change very slightly in the doping range 
[0.05-0.4]. When x exceeds 0.4, though, region II be- 
comes increasingly thinner and localized around the V 2 
axis. The only role of Vq on the d-wave state is to set its 
domain of stability. In particular, the magnitude of 
does not depend on it. 

The size of the gaps as a function of hole concentration 
for pairs (Vb, V 2 ) inside regions II, IV and V is plotted in 
Figure 2. The three figures taken as a whole show that 
the d state is more stable only close to half-filling and for 
large enough values of the ratio V2/V0, while the s states 
exist in the whole range of hole dopings. The shape of 
the curves for is independent of Vq; V 2 only changes 
their height, so that one can write A^ = f(V 2 )g(x). They 
always reach a peak at x — 0.2 and vanish at x±0.7. The 
form of the curves for the s-state gaps does depend on 
both Vq and V 2 , on the other hand, but its two maxi- 
mums always occur at x = 0.2 and 0.8. We find that 
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FIG. 2. Gaps (measured in meV) for (a) Vo = —250, 
V 2 = 150 meV; (b) Vo = 50, V 2 = 150 meV; and (c) V = 250, 
V 2 = 50 meV. 

the gaps of the t\ — t 2 model have peaks at the same 
doping concentrations. We therefore conclude that their 
existence and position is due purely to filling effects and 
not to the peculiarities of the band structure. The arrow 
in Fig. 2(c) marks the boundary of stability x c between 
the d and the s states, the free energy of the d state being 
lower for x < x c . 

There are compelling reasons to believe that the renor- 
malized on-site interaction of this low energy model 
should still be strongly repulsive, so we focus our study 
on the negative Vq quadrant now. A first point to notice 
is that the d-wave superconducting solution disappears 
for values of the on-site effective repulsion of the order of 
two times Vz- This fact will set some constraints on the 
possible parent Hamiltonian of this model if it survives 
the inclusion of fluctuations over the Mean Field (MF) 
solution. Moreover, we estimate that the strength of the 
effective nearest neighbor attraction needed to obtain the 
experimental value of the gap for optimally doped Bi2212 
(30 meV) is of about 140 meV, while Vo can range from 
-320 to meV. Such values, compared to a band width 
of 1.1 eV, place the resulting model within the weak cou- 
pling regime. This result leads to question the quanti- 
tative-Accuracy of the BCS gap and of the whole MF- 
RPAlj scheme at intermediate couplings even at T = 0. 
What actually happens is that such theories provide fine 
estimates for the gap and other physical magnitudes for 
large values of the coupling constants while for interme- 
diate couplings and low dimensions, the results are much 
poorer. A comparison of MF theory and Bethe ansatz re- 
sults for the Id Hubbard model at half filling and T = 
can be found in Ref. where it is shown that the size 
of the gap is grossly overestimated by the MF solution at 
intermediate couplingrj We think that the gap should 
be estimated beyond MF-RPA theory in order to make 



TABLE I. Changes in the gaps due to orthorrombicity for 
the case Vo = —250 meV, V2 = 160 meV and x = 0.2. Gaps 
are measured in meV. 
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direct comparisons with experiments. Such a theory will 
require much larger values of Vo an d Vi to give a gap of 
30 meV, thereby placing the coupling constants in the 
intermediate coupling regime. 

The ratio 2Ad(kF, ma x,T = Q)/KbT<i as a function of 
doping is non-universal and larger than the value ob- 
tained for a parabolic band (4.14). There is a plateau, 
which ranges from x = 0.05 to 0.4, where the ratio is 
almost constant and equal to 4.35, while it grows for 
smaller or larger doping, even reaching a value of 7 for 
x = 0.65. 

The present model provides a universal value xq — 0.2 
for optimal doping, defined as the hole concentration 
which yields the maximum critical temperature. We 
should caution however that experimentally, Ad{T = 0) 
and the temperature T* at which the pseudogap begins to 
open up are monotonically decreasing functions of dop- 
ing and do not peak at xq, while the experimental. X 1 ^ h as 
a maximum at about x — 0.15.E3 RPA theoriesE2lE3 of 
s-wave superconductors, which permit to study the inter- 
mediate and strong coupling regimes find that electrons 
begin to pair up at a temperature T* which is higher 
than the critical temperature T c where pairs condense. 

We study now how orthorrombicity effects can induce 
mixing of d and s states as might happen in YBCO, a 
cuprate where the 2. and b crystal axes become substan- 
tially incquivalcnt.H In order to do so, we redefine the 
hopping integrals ti tXtV = i;(l ± a) for i = 1, 3, 4 so that 
the probability amplitude for electrons to leap between 
two sites depends on the spatial direction of the hop. We 
also suppose that the orthorrombicity leads to a different 
nearest neighbors interaction in the x and y directions, 
V~2,x,y — ± j8). We find that a finite j3 doesn't mod- 
ify much the d-wave gap, but indeed induces a finite s 
component (see table I) . Orthorrombicity in the hopping 
integrals produces appreciable changes in A^, but leads 
to smaller modifications in the s-state gaps. Notice finally 
how a sign change in a and/or [3 inverts the sign of the 
s components. This effect, which was already deduced, 
phenomenologically using Ginsburg-Landau theoryj£3'Ej 
implies that the s-wave component of the gap has oppo- 
site sign in different twin domains of twinned samples of 
YBCO. This is also the explanation given for the Fraun- 
hofer pattern observed in a recent c-axis Josephson tun- 



3 



ACKNOWLEDGMENTS 



> 



o 



5 
4 
3 
2 
1 

-1 
-2 
-3 



* 

» 


: n 8 | (b) : 
o ■ 


8 

S 


■ 


• x=0.02- 


• x=0.02 s 


- x=0.20 , 


■ x=0.20 


♦ x=0.30 


♦ x=0.30 9 


* x=0.50 


* x=0.40 


x 0.60 ! 


x=0.50 



4 6 



4 6 8 



1/V (eV ) 1/V 2 (eV ) 



FIG. 3. (a) log As versus 1/Vo for Vi = and several x. 
(b) log A d versus I/V2 for V = 0. 

neling experiment, where the junction had been grown 
across a single twin boundary.Ej 

One may compare the results obtained within this 
tight-binding model with the predictions of a much sim- 
pler isotropic Hamiltonian, with a linearized dispersion 
relation and pairing interaction of the form 



6 = 
V{k±k') 



2m 



— fl ~ T%VF(k — kp) 

V Q L - V 2 L cos(2(9) cos(26>') 



(5) 



In this case, the gap function A L = Afp-W A d cos(26*) 
obeys the following saddle point equationstS 



g L v , 



(16 



de — , ; ... — — , (6) 



y/e 2 + \A L (0)\< 



where F S L = 1, F£ = 008(20), and g L is the two- 
dimensional density of states. Now, expanding Eq. (Q) 
up to order (£;£) 4 and comparing with Eq. do), we obtain 
the following relations: g L V$ — gVo, g L V} = j 2 gV 2 , 
= A s and A d = j Ad, where g is an effective tight- 
binding density of states and 7 = -(/cF£) 2 /2 + (fcF0 4 /24. 
These two parameters, together with the cut-off E4, can 
be used to fit A s ,<j with A^ d along the Vo and V 2 axes, 
respectively. As shown in Fig. 3, we indeed find linear 
behavior for log A as a function of 1/V, from which we 
extract values of /cfC/ 71 " ranging from 1.21 to 1.27, and g 
from 1.1 to 1.6 eV~ l . 

In conclusion, we have proposed and studied an effec- 
tive tight-binding model for the superconducting state 
of the cuprates, guided mostly by results of ARPES ex- 
periments. In particular, we have delineated its weak- 
coupling phase diagrams for different hole doping con- 
centrations close to the experimental optimum doping. 
We have been able to establish the region of parameters 
which is likely to be relevant for the cuprates in absolute 
energy units (meV). We have finally mapped the tight- 
binding model onto a much simpler linearized model and 
obtained similar phase diagrams. 
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